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Type-I matrices were introduced recently as finite dimensional prototypes of quantum integrable 
systems. These matrices are linearly dependent on an "interaction" type parameter, and possess 
interesting properties such as commuting partner matrices and generically violate the von Neu- 
mann Wigner non crossing rule. The important role of Pliicker relations in this construction is 
noted. Type-I matrices are given a transparent formulation in terms of Fermi or Bose type particle 
f"^ | operators- they represent a Quantum glass model with either Fermi or Bose statistics, with several 

, free parameters that may be chosen at will. 
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INTRODUCTION 



In Ref [l[ (referred to as I below), the author initiated a study of quantum integrable systems in finite dimensions, 
within the context of parameter dependent commuting matrices. This in turn was motivated by several examples 
of specific integrable models, such as the Hubbard model0-0| and the Heisenberg model [jj. In these examples, 
one studies the realizations of the general model in Fock space for particular sectors of quantum numbers, such as 
momentum, parity, total spin and number. This lead one to real symmetric matrices in various dimensions N > 2. 
These have the remarkable feature that the von Neumann Wigner non crossing rule@ is violated. One ends up with 
several level crossings that are conventionally termed "accidental" . This terminology is rather avoidable, since there 
is a also a belief that there is nothing accidental in having such level crossings; the existence of several dynamical 
conservation laws (dependent on coupling constants) are believed to be causally implicated. Further, the statistics of 
energy levels of these integrable models are also known to be close to Poisson statistics, and hence consistent with 
£H . the absence of level repulsion that generic systems are known to possess 0, Q. 

While the current general programme for the study of quantum integrable systems focus on properties such as 
factorizable S matrices, or the Yang Baxter relation, the approach of (I) Ref [l[ gets to the core of the issue of the 
matrix realizations of these models. As such, it is "blind" to the specific physical details of the models. The main 
results of (I) Q is the identification of a class of matrices, termed Type-I matrices, discussed in detail below. Here 
the core property of multiple parameter dependent conservation laws is made explicit, and one has an algorithm for 
generating such matrices as well as a count of the number of such matrices. 

In an impressive work, Ref.@ Owusu, Wagh and Yuzbashyan (OWY) have built on this initial advance, and 
£f~) • produced several further results. A fundamental advance is the introduction of a basis of matrices, in terms of which 
the matrices of Ref. [l[ can be expanded. OWY further show a link with an integrable model due to Michel Gaudin [Io| , 
that is currently enjoying popularity in the context of superconductivity of finite systems jioj. OWY also throw light 
on the "mechanism" of the the level crossing, and give explicit formulas for the number of level crossings one finds in 
Type-I matrices. 

The objective of this paper is multifold. Firstly, a fundamental constraint equation in Ref.[l[ for constructing 
Type-I matrices is shown to be related to the so called Pliicker relations of Pfaffians. Since Pfaffians are basic to 
anticommuting objects such as Majorana Fermions, one sees that Fermi statistics enters this program of describing 
integrable systems in a fundamental and unexpected fashion. From this analysis, the parameterization of the solutions 
of Type-I matrices by OWY in Ref.jij], arises as an elegant consequence, and the entire construction becomes more 
transparent. 

Secondly, I show that the link with the Gaudin type model[Io| is made more naturally with Fermi (or canonical 
Bose) statistics. The connection made by OWY with the Gaudin model assumed hard core Bose statistics for the 
particles, and is confined to the sector of one spin wave, i.e. is confined to a specific sub manifold. The basic quantum 
operator underlying this class of problems is the permutation operator that has several possible realizations, leading 
to distinct models. The permutation operator has a Fermi representation: this is shown to be more natural than 
the (hard core) Bosonic one used in Ref.[9|. Once the commutation relations of a set of matrices is established, we 
can elevate these to operator relations with either Fermi or Bose statistics (see Eq. (fT9|) below) and thus also have a 
Bosonic representation of these. 

With either Fermi, or with soft core (i.e. canonical) Bosons, we construct a Quantum glass model below, i.e. a 
particular type of Anderson model for disordered carriers. This model is akin to Gaudin's model with hard core spins, 
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and depends on several parameters that may be chosen as one wishes, and has commuting partners in all particle 
sectors. These commuting partners may be thought of as local charges that are broadened out in a specific way. 

Interestingly the Phicker relations arise in other aspects of integrable systems as well, these are central objects in 
Sato's work on classical solitonic theories (i.e. classically integrable systems) [ll, 12 1, where the so called r functions 
satisfy these relations. For quantum integrable models, a connection has been shown to exist between the transfer 
matrices of and the bilinear identities of the r functions 13[ satisfying Pliicker relations. 

A few remarks are useful to put the current work and the related Refs. [l], Q , within the context of matrix theory as 
used in Quantum theory. In order to keep things simple, let us specialize to finite dimensions |2l|. Quantum observ- 
ables lead to Hermitean, or real symmetric matrices, and simultaneous measurability of observable pairs translates to 
the theorem that commuting matrices of the above type are simultaneously diagonalizable. One simple result about 
two such commuting matrices a and b, is that one of them is expressible as a power series in the other p2j . The 
current series of works differ from these in that the focus is on matrices that depend in a simple way (linearly, or 
possibly a polynomial of low degree) on a parameter, and one insists upon the commutation property for all values of 
the parameter. This problem is natural in the context of examples in Quantum theory, such as the hydrogen atom, 
where the Laplace Runge Lenz vector depends linearly on the squared electric charge. It is also true in the structure 
of the higher conservation laws in models such as the Heisenberg [H[ and Hubbard models 0, Q|. The notable results 
in the works Refs.[l|, |9] follow from the detailed study of the simple parameter dependence of the commuting pair. 



SUMMARY OF (I) AND THE INTRODUCTION OF A BASIS OF COMMUTING OPERATORS 

In Ref [H, we introduced a family of real symmetric matrices in N dimensions depending linearly on a parameter x. 
These were introduced as purely algebraic prototypes of integrable systems in finite dimensions and termed as Type-I 
matrices. They are efficiently represented as: 

a = a + xA, with A = A<j + [a, S], (1) 

with tw o g eneric diagonal matrices a and A^ having unequal entries, i.e. a = DiagonalMatrix{u\, . . . ,itjv} with 
Ui 7^ Uj [23j| . Ad = DiagonalMatrix{Ai, . . . , An} with A{ ^ Aj, and a real antisymmetric matrix Si.j — —Sj t i- One 
may think of a and A as the kinetic and potential energy matrices, and the parameter x as a perturbation parameter 
in typical quantum systems. In this notation, the role of the antisymmetric matrix S%j is made explicit. There is 
no loss of generality since if we are given the matrix a(x) in an arbitrary basis as the sum of two non commuting 
matrices, we can convert it to this form by performing an orthogonal transformation that diagonalizes the matrix 
a(x = 0). 

In this way, we model integrable systems, without reference to their explicit origin in the physical world, as parameter 
dependent matrices. This construction is inspired by the standard examples of the Hubbard and Heisenberg models. 
In these models, finite dimensional matrices of the above type emerge on restricting the state space to various sectors 
of usual (parameter independent) conservation laws such as particle number, parity, spin and total momentum. 

Since integrable systems are known to possess several parameter dependent (i.e. dynamical) conservation laws, one 
wants to know if other matrices depending on x, possibly linearly, can be found. It was indeed shown that under 
certain conditions on S, summarized below, such commuting partners (3{x) can be found, i.e. [a,/3] = 0. The form of 
the dynamical conservation laws fi{x) was shown to be very similar as that of a(x): 

/3=b + xB, with B = B d + [b,S], (2) 

where b and B^ are diagonal matrices. 

In (I), it was shown that the number of matrices of type-I in N dimensions is V a = (3N — 1), and for a given matrix 
a(x) out of this set, there are a further V& = N + 1 matrices of the type (3(x). 

The crucial condition on the antisymmetric matrix S was written in (I), in terms of its inverse elements Ri.j = 

4>(i,j,kd) = Ri,jRk,l — Ri,kRjA + Ri,lRj,k = 0. (3) 

These equations are extensively discussed in mathematics literature as the Pliicker relations (l4| . and their analysis 
is presented later. In our original work (I), we noted that these are greatly overdctcrmincd equations, since there are 
N Ci quartets of indices and equations, but only C2 matrix elements Rij to be determined. In (I) it was shown, by 
using a consistency condition involving 5 indices(Eq(I-15)), that this set has Vr = 2N — 3 free parameters and hence 
independent solutions. For example one may choose at will the parameters R\.j\ 2 < j < N and i?2,fc! 3 < k < N, 
and the remaining Ri m are determined in terms of these with no conflicts. 
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In Ref.@ Owusu, Wagh and Yuzbashyan (OWY) have shown that it is more efficient to introduce a basis of 
commuting operators in terms of which both the matrices a Eq. ([1]) and /3 Eq. ([2]) can be expanded linearly. Here a 
commuting basis guarantees the commutation of the matrices a and /3. The basis of commuting operators {Z(r)}, 
with 1 < r < N may be written in terms of the Dirac projection operators as 

Z(r) = |r)<r| + xj^i P.(r)\s)(s\ + S r>s (\r){s\ + \s){r\) } . (4) 

S 

Here S r , s will be seen below to be the same elements as in Eq. ([1]). The following commutator vanishes: 

[Z(r),Z(s)]=0, (5) 
provided the matrix elements Pi{j) satisfy the conditions 

5>.(»0 - 0, (6) 

r 

Z>.(»o = °. ( ? ) 



A(i,j;l) = Pi (l)- Pj (l). (8) 



The constraint Eq. ([7]) guarantees that the trace of the operators vanishes, it is not necessary for the commutation of 
the Z's but is a convenient condition. 

Now Eq. ([8]) can be rearranged in a way that eliminates the pi(j) variables as a four index identity (where "1" is a 
spectator index) 

A(*,j;0 + A(j,fc;0 + A(*,*;0 = 0, (9) 

and in this form it is identical to Eq. (1-10) , and by using the inverse matrix elements Ri.j = jh-, it becomes precisely 
Eq. ([3]) above. The point of this construction is that we can now take sums of the basis operators in Eq. (fTJ) 
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and in this way recover the matrices found in (I) . By allowing some of the cij to be pairwise equal, OWY obtain a 
somewhat greater freedom than that in (I), where all the Uj were chosen to be distinct in order to obtain generic 
matrices. 

The number of independent parameters in the S's or equivalently in R's is N(N — l)/2, and as mentioned above, 
we showed in (I) that the constraints in Eq. ([3]) are mutually consistent, giving 2N — 3 free parameters in S. OWY 
parameterize the solutions of Eq. by a neat ansatz, namely 

Rl , 3 = iii^M. (10) 

These are obtained in turn by performing a local gauge transformation Eq. (|11[) . on a particular solution Ri j = 
{si — Ej) noted in (I). This transformation consists of multiplying each matrix element by an arbitrary j dependent 
factor, 

Rij -* V(7i7j)-Rt,i) <Pijkl <t>ijkl /{liljlkll), (11) 

and is clearly a way of generating further solutions from a given one. 

Indeed Eq. (flUl) has the correct number of parameters (27V — 3). To see this, we start with the N e's, and the N 
7's giving us 2N parameters. As discussed more fully below in Eq. (fT7|) . we subtract 3 parameters from 2N, since 
we can shift all e's by a single constant and further scale all the e's and all the 7's by two j independent constants. 
Using Eq. ([8j), we may then infer the Pj(i) from this parametrization of R^ — 1/&j, and find 

7 2 

Piti) = — and 

Pi(i) = 7, 2 E^— ( 12 ) 
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We will use this convenient parameterization in the rest of this work. For completeness, we note that the parameter- 
ization of Eq. ([1]) i.e. a — J2 u jZ(j) in (I) translates to the new variables as follows: 



Sij 






£i - e s 






= const — 












if 

£ i ~~ £ j 


k^i 



(13) 



with Yij from Eq. (1-5), so that n(i;jk) in Eq. (1-19) vanishes identically. 

PLUCKER RELATIONS AND THE PARAMETRIZATION OF THE ANTISYMMETRIC S MATRIX 

Eq. ((3]) were recognized belatedly by the author, as Pliicker relations of mathematical literature. Since these are 
central to the construction of this class of matrices we take a closer look at the 2N — 3 solutions that were found in (I) . 
We explore the structure of the relations by using a more rigorous technique next, and see that the ansatz of OWY 
follows from the analysis as the unique solution. We show that these relations involve the so called Pliicker relations 
for Grassman variables, and hence presage the final form our presentation that involves Fermions in a fundamental 
way. 

We begin by noting that EqQ involves 4>Ujki)i which is a Pfaffian of a real 4x4 skew symmetric matrix Rij. 
The vanishing of <f>Ujki) is a standard example of a Pliicker relationjl4j. The totality of these equations is expressed 
elegantly using exterior forms. Let us define a N dimensional real vector space spanned by unit vectors ej and define 
an antisymmetric wedge products e, A ej. These provide a basis for the linear vector space W(2)[14|. In this space, 
we define for a skew symmetric Rt.j a "two form" : 

1Z = Rij ei A ej. 

i<j 

It is now easy to see that 

TZATZ = 2 ^2 <i>i,j,k,l &i A ej A e k A e u 

i<j<k<l 

and hence we recognize that the totality of relations in Eq. ((3]) are precisely equivalent to finding solutions of 

KAK = 0. (14) 

This condition defines [3] the "decomposability" of the two- form 1Z. This problem can be resolved by noting that 
every skew symmetric matrix can be expressed in its real normal form involving 2n orthonormal vectors a a and 
b a with 1 < a < n, satisfying the conditions J^j Ri.j a< j — A Q 6f and ^2,jRi,jVj = —\ a af. We may term these as 
the pseudo eigenvectors and pseudo eigenvalues, since the Hermitean matrix i R has real eigenvalues ±A" and real 
eigenfunctions -^(dj ± an d n < N/2 is the number of non zero eigenvalues of i R. The normal form is expressed 
as 

/•',.. = E A " ":'''.'/• 

a — l,n 

With this decomposition, and with v a = V^"J2j a< j e 3 ano - u;Q = "V^" J2j &j e ji we can rewrite the relation 

K= \ V v a Aw a . 

a=l,n 

We thus see that Eq (|14[) is possible if and only if the number of vectors n = 1, i.e. there is only one pseudo eigenvector 
of R. This is known as the condition of decomposability, and provides us with a neat representation EqQ with a 
single eigenvalue A and a pair of orthonormal vectors Xj y, j 



Ri,j — ^{XiUj — ViXj)- 



(15) 
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Using the local gauge invariance in Eq ([TTj) . we can drop the condition of orthonormality of Xj and yj in Eq (fT5]) . as 
far as generating solutions to the original problem EqQ is concerned. We may also absorb the A" factor into the 
vectors, and it appears that we have 2N independent real parameters in the solution of EqQ. However, we observe 
that there is a redundancy in this counting, the vectors Xj and yj can be changed without changing Ri j if use three 
linear transformations with arbitrary parameters p, q, r as 

{xj,yj) -> (xj +pyj,yj) 
{xj,yj) -> (xj,yj+q xj) 

(xj,yj) -> (r Xj, - yj). (16) 

We thus see that the total number of real parameters available is exactly 2N — 3 as known already from (I). One 
convenient set of 2N — 3 variables was given as R\ j with 2 < j < N, and i?2 ,j with 3 < j < N, in terms of the 
Xj,yj we may e.g. set x\ = l,y\ — l,x% = 1 and determine the remaining 2N — 3 variables from the -Ri.j's. The 
parameterization EqlfTTjj) of OWY can be obtained from Eq ([T5|) by setting Ej = ^ and 7j = ^- , and the symmetries 
of Eq tfTB")) are transformed into 

(£j>7j) ~> + P,7j) 

(ej,jj) -> (r 2 ej ,r 7j ). (17) 

We may again reduce the apparent 2N parameters by 3 using these relations, it amounts to choosing three parameters, 
say (£1,71,72) arbitrarily as e.g. (1, 1, 1) and the rest are fixed using the inverse of Eq (fPD"| . 



FERMIONIC REPRESENTATION OF COMMUTING OPERATORS 

We next show that the matrices Z(r) in Eq. (U) lead to a neat Fermionic representation, which may be thought of 
as a model for a Fermi glass with localized states. Let us define a Fermionic set of operators flj,a] and rij = aja^, 
obeying the standard anticommutation relations 

{ai,aj} = 5ij, (18) 

with 1 < i < N. It is elementary to see that two commuting matrices [P, Q] = lead to a commuting set of Fermionic 
operators (e.g. see [HI), i.e. 

Pij a i a j: Qij^j] = ^[P^h-majam = 0, (19) 
ij ij Im 

where [P, Q] is the matrix commutator of the two matrices and Q,j . Thus we obtain a set of TV Fermionic operators 

Z(r) = n r + [ p s {r)n s + S r>s (aja s + a\a r ) ] . (20) 

We see that these inherit the commutation property [Z(r), Z(s)] = from Eq. Using the parametrization Eq. (flU)) 
and Eq. (|12p , we write the basis set of commuting operators as 

- 1 

Z(r) = n r + x*^2 — - — [7r7s (ata s + a\a r ) - 7^ n s - 7^ n r ] , (21) 

s e r e s 

where the prime indicates s 7^ r. 

Readers wishing to skip the earlier discussions, can directly verify that the commutator [Z(r),Z(s)] vanishes, for 
arbitrary values of the given parameters by a straightforward calculation. 

We also remark that the choice of the statistics of the canonical operators a,j is not the only one possible. The 
entire argument of this section can be repeated if we use canonical Bosonic operators instead, i.e. aj — > bj where 
[bi, bj] = 6ij. Thus one can equally well consider a Bosonic glass rather than a Fermi glass model here. 

Finally we note that the single particle sectors of the Bosonic, Fermionic and hard core Bosonic models are all 
identical and correspond to Type-I matrices. For higher numbers of particles, these correspond to other classes of 
matrices depend on the statistics chosen, e.g. these are Kronecker products of Type-I matrices in the case of canonical 
Fermions and Bosons. 
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Mapping to The Gaudin Model 



The mapping discussed by OWY views Eq. (|4]) as the S^ otal = N/2 — 1 subspace representation of the Gaudin 
Hamiltonian 3| 

: E 7^-^-4 (22) 



7Gaudin 



= Sf 



Si - £j 



This model was first written down by Gaudin[10]. Gaudin actually wrote it without the first term Sf, this was 



supplied later by Sklyanin 16[ from twisting the boundary conditions. To be exact Eq. fU) has an extra factor of 
7i7j that OWY argue can be incorporated into the equations, and also their magnetic field term B is cx l/x. The 
Gaudin model is currently very popular for describing the dynamics of Cooper pairs within the BCS theory for finite 



systems [17|,ll8|. Each spin flip represents a Cooper pair, from the Anderson mapping of the BCS theory to spin waves. 
Thus = Ckiic^k^ and the label i is actually a momentum space label. 

The point about the Fermionic representation Eq. (|2ip of Eq. Q is that it is true for all numbers of Fermions, and 
not restricted to a single particle sector. In this sense, the present Fermionic representation is much more powerful, 
and further the factors ji do not need any special treatment, they are automatically treated in the commutation 
relations. Thus Eq. (|4j are embedded without any further qualifications in the operator equations Eq. (|2"Tj) . We see 
below that this representation enables us to find applications of this model for Fermions in a disordered potential, i.e. 
the Fermi glass problem. The same statement is also true if we use canonical Bosonic operators instead of Fermions, 
as mentioned above. However the Gaudin model is expressed in terms of hard core Bosons, and the magnitude of the 



spin is related to the jj making the scheme somewhat cumbersome [24 1 



DIAGONALIZING THE FERMI HAMILTONIAN AND THE DENSITY OF STATES. 

We now turn to a study of Eq. (|21[) and a related Hamiltonian obtained by summing 

Hr = ^ e r Z{r) = ^ n r£r + x ^ lilja\a _ x N X! 7 J ' ( 23 ) 

r r ij j 

where N — ^2 r n r is the number operator. For Fermions or canonical Bosons, this Hamiltonian is the analog of the 
so called Richardson [l9j| Hamiltonian in the theory of nuclear matter (the 7 factors do not usually arise in the latter). 
The Richardson model and also the related BCS[l7| problem for finite systems [l^, are expressed in terms of hard 
core bosons (i.e. spin half objects) representing Cooper pairs S~ = c^C—j^f . These are in turn, obtained by taking 
sums over the Gaudin Zj }audm operators of Eq. (|22p . The Hamiltonian Eq. (|2"5)l is considerably simpler to solve for 
a general population of particles than the corresponding problem for hard core bosons, and is akin to a free gas of 
particles in a suitable one body potential. Clearly our Hamiltonian Eq. (|2"5|) commutes with each of the Z(r), and 
plays a central role in the Fermi glass interpretation. 

If we view the labels i,j as wave vector indices, then H describes a band model with an arbitrary dispersion £j. It 
is subject to a potential that scatters from every wave vector to each of the others, with a potential matrix element 
X7i7j. Since the "fi are arbitrary, they may be chosen at random. We thus realize a band model with a separable 
random scattering potential. If on the other hand, we view i as site labels in a tight binding model, the energies e, 
may be chosen at random, and the kinetic energy hops between every pair of sites- i.e. realizing an infinite ranged 
random Fermionic Anderson model. 

We now turn to the task of diagonalizing the Hamiltonian and all the Z(r) by an ortho gon al transformation. This 
transformation for the single particle sector is essentially identical to the one in Richardson [19] , and many subsequent 
works, and hence we will be brief. Define new canonical Fermion set 

4 = J2 Q ^ 4 with = 5 «> ( 24 ) 
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through an orthogonal transformation generated by a real orthogonal matrix Q such that Q T ■ Q = 1, and 

7i 



E 



Wj(x) 



(25) 
(26) 



- = £— — • ( 27 ) 

X U) m (%) - £j 

J 

Here uj m (x) in Eq. (f27|) are the iV eigenvalues of the transformed Hamiltonian Eq. (|23|) . we write the argument (x) 
to emphasize that these depend parametrically on x. A short further calculation gives 

H R = "mW 4^ m - x N 7j 2 , (28) 



J 



m =x ^^-^— d i dm . (29) 



As x — > 0^, one sees that ojj — > £j±0, i.e the w's are pinned to the e's. The eigenvalues u) m {x) are in 1-1 correspondence 
and evolve out of the numbers e m smoothly as x increases from zero. Thus the eigenvalues ujj of the Richardson 
Hamiltonian Eq. (|23p interlace the numbers Ej, with one extremal eigenvalue that grows linearly with x. For a: 3> 
(x <C 0), the extremal eigenvalue ujn ^> £n (^i <C £i). The density of states of Wj has a width that remains fixed 
with x if we ignore the exceptional extremal case. It is easy to see that the uj m do not cross each other as x varies, 
and they do satisfy the von Neumann Wigner non crossing rule. The conserved quantities Zj may be visualized as 
evolving continuously from the occupation numbers rij as x increases from zero. 
We could more generally consider the two operators formed from the sums 

a{x) = J2u r Z{r), and /3(x) = ^ v r Z(r), (30) 

r r 

with arbitrary u r and u r , and see immediately that these are the Fermi space representations of the operators 
introduced in Eq. ([TJand Eq. ([2]): 

a(x) = 2J u r n r + - ^ — [iris {a\a s + a\a r ) - 7^ n s - 7^ n r ] 

r r,s 
I 

K x ) = E Vr 7lr + 9 E ' _ [iris (ata s + ala r ) - 7^ n s - 7^ n r ] . (31) 

These commute mutually for any choice of the parameters, including x, and also with the constants of motion Zj in 
Eq. (|2"Tj) . and on diagonalization become 



a(x) = y^q m (x) d^dm, 



9 

A comment on the conservation laws Eq. (|2ip and their relationship with the "Hamiltonian" Eq. (|3~Tj) is useful here. 
At x = the existence of N constants of motion of the Hamiltonian is obvious since the Z' s are just the number 
operators of the Fermions. When we perturb the Hamiltonian from this "free case" by adding any term proportional 
to x, we can always fix the conservation law to be valid to 0(x), but generally the terms do not commute to 0(x 2 ). 
This is familiar in the theory of integrability violating perturbations to integrable systems, as in the Kolmogorov- 
Arnold-Moser theory [20|: the conservation laws analogous to Eq. (|21j) can be rescued to linear order in the new 
perturbations, but not to higher orders. The speciality of the specific perturbation in Eq. pip is that there are no 
corrections to 0(x 2 ) and the the conservation law Eq. (|2"Tj) commute exactly. 
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The inevitability of level crossings for Type-I matrices was noted empirically in (I) , on the basis of several examples 
that were studied. However the theoretical explanation awaited the work of OWY, who showed that for a generic 
choice of u r , the eigenvalues of a(x), i.e. a m {x) in Eq. have atleast one and at most N ~ x C<z level crossings as x 
varies over its range. These eigenvalues thus defy the von Neumann Wigner non crossing rule, unlike the eigenvalues 
of the Richardson Hamiltonian w m (i), which do obey the rule. The one exceptional case is Uj — £j when the a(x) 
reduces to the Richardson hamiltonian Eq. (j2"5)l 25| . One may understand the violations of the non crossing rule by 
thinking of the eigenvalues of d(x) as smeared versions of ui m {x), and thereby less sharply governed by the rule. This 
is illustrated in Fig. 1., where we plot the energy levels for N = 5, and show that while u) m (x) of Eq. (|27[) avoid level 
crossings, the derived eigenvalues a m (x) from Eq. ([30l l32|) do display level crossings. In this sense, there is a hidden 
generic model Hr satisfying the non crossing rule, behind the violations of the same in the constructed matrices a{x). 



FIG. 1: Left panel shows the eigenvalues u m {x) from Eq. (127[) for the case of N = 5, where one observes narrowly avoided 
level crossings involving the top three levels. Right panel shows the effect of mixing levels through Eq. (|30l 132 with Uj — Sj 
chosen randomly with a small scale of variation. We see that the eigenvalues a m (x) obtained from Eq. (|3UI32p . cross each other 
profusely, thereby violating the Wigner von Neumann non crossing rule. 

Finally, we note that the Hamiltonians Eq. (f5T|) with Fermionic ctj (Bosonic bj) can be viewed as representing a 
class of localized states in the Fermi (Bose) glass problem of disordered non interacting quantum particles. At x = 0, 
the model consists of localized states with energies u r , and clearly has N conservation laws Z(r) as in Eq. (|21|). 
corresponding to the occupation numbers of the different sites. As x varies from zero, the particles hop around as 
dictated by the Hamiltonian, but with generalized conserved occupancies at all sites given by Eq. (|2~T1) . These are 
therefore localized to all orders in the perturbation x, despite hoppings that carry them far away. We can easily see 
that the energy level statistics of these systems follow the Poisson distribution for small separations, due to the level 
crossings that occur in these Hamiltonians. The absence of level repulsion what one expects from localized states in 
the Anderson model on general grounds. 
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